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ABSTRACT 

We investigate the pulsation driving mechanism responsible for the long-period photometric varia- 
tions observed in LS IV-14°116, a subdwarf B star showing a He-enriched atmospheric composition. 
To this end, we perform detailed nonadiabatic pulsation computations over fully evolutionary post 
He-core-flash stellar structure models, appropriate for hot subdwarf stars at evolutionary phases pre- 
vious to the He-core burning stage. We found that the variability of LS IV- 14° 116 can be attributed to 
nonradial g-mode pulsations excited by the e-mechanism acting in the He-burning shells that appear 
before the star settles on the He-core burning stage. Even more interestingly, our results show that 
LS IV-14°116 could be the first known pulsating star in which the e-mechanism of mode excitation is 
operating. Last but not least, we find that the period range of destabilized modes is sensitive to the 
exact location of the burning shell, something that might help to distinguish between the different 
evolutionary scenarios proposed for the formation of this star. 

Subject headings: Stars: oscillations — evolution — subdwarfs — individual (LS IV-14°116) 



I. INTRODUCTION 

About 5% of all hot subdwarf stars (sdB, sdO) 
show He-enriched surface abundances (He-sdB, He-sdO). 
While most normal H-rich subdwarfs are supposed to be 
low mass core He-burning stars with atmospheres domi- 
nated by H due to the action of gravitational settling, the 
evolutionary status of the He-rich subclass is less clear. 
He-rich subdwarfs have been sugg ested to be the result 
of either He white dwarf mergers JSaio fc Je ffeiTl[2000l) 
or late helium core flashes (j Brown et al. 2001). Regard- 
less the particular scenario proposed for their formation, 
it is accepted that some of these stars are still contract- 
ing towards the He-core burning phase (EHB), as other- 
wise gravitational settling of the remaining H would have 
turned the star into a H-rich sdB star. 

Two families of pulsators have been discovered 
within the H-rich sdB st ars: the rapid pulsators (sd- 
BVr; Kilk ennv et al.|[201 0l) discovered by Kilkenn v et al 



(|l9 97), and the slo w pulsators ( s dBVs ; Kilken nv et al 
[201 0() discovered bv iGreen etall (|2003D . While sdBVr 
stars show short pulsation periods 80 — 400 s) as- 
cribed to radial modes and non-radial p- modes, sdBVs 
pulsations (with periods ^ 2500 — 7000 s) are associ- 
ated to non-radial long period g- modes. Pulsations in 
both groups of variable stars have been explained by 
the action of the K-mechanism due to the partial ion- 
ization of iron group elements in the outer layers, where 
these elem ents are enhanced by t he action of radiativ e 
levitation (jCharpinet et al.l 119971: [Fontaine et alll2003[ ). 
While many sdBV stars are known, only one He-sdB 
star, LS IV-14°116 , has been f ound to be a pulsator 
([Ahmad et al.|[200l [Ahmad fc J efl^crv 2005). In fact, LS 
IV-14°116 is a very intriguing object: while its atmo- 
spheric parameters (Tgff and q) place it well within the 
sdBVr instability region ([Naslim et al.|[20Tol : IGreen et al.l 
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1201 If ), it dis plays periods typi cal of the sdBVs family 
of pulsators jAhmad fc Jefltervll2005l: IGreen et aL| [2011). 
From a spectroscopical perspective, LS IV-14°116 is 
also an intriguing object, showing a mild He-enrichment 
(riHc — 0.16) and very high excesses of z irconium, yt- 
trium and strontium in its atmosphere (jNaslim et al.l 
[20Tl[) . In this sense, both the He-enrichment and the 
heavy metal excesses might be related to the effects of 
ongoing diffusion before reaching diffusive equilibrium, 
thus placing the star in the pre-EHB phase. Up to now, 
the driving mechanism behind the long period pulsations 
of LS IV-14°116, as well as the absence of short period 
pulsations, remains a mistery and the st ar poses a chal- 
lenge to the theory of stellar pulsations ([Fontaine et al.l 
120081) . In this connection, the recent confirmation of both 
the multiperio dic variability and the Tes — g values for 
LS IV-14°116 (IGreen et al.ll201lD strongly increase the 
enigma. 

While most pulsations exhibited by pulsating stars (in- 
cluded those of sdBVs and sdBVr stars) are self ex- 
cited through the classical K-mechanism, no star has 
been discovered so far to be excited by the e-mechanism. 
The only pos sible exception are some oscillations in 
5 Scuti stars ( Lenain et all 120061) and PG 1159 stars 
([Kawaler et al.lil986HG6rsico et al ll2009[ ). but neither of 
these have been confirmed. In the e-mechanism, the driv- 
ing is due to the strong dependence of nuclear burning 
rates on temperature. In the layers where nuclear reac- 
tions take place, thermal energy is gained at compres- 
sion by the enhancement of nuclear energy liberation, 
while the opposite happens during the expansion phase 
(lUnno et al.lfl98l . 

In the present letter we show that long period g-mode 
pulsations, like those observed in LS IV-14°116, can be 
understood as driven by the e-mechanism acting during 
the unstable shell burning events that take place before 
the star settles on the He-core burning phase. In addi- 
tion, our scenario might explain the intriguing absence of 
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TABLE 1 

TiPICAL TIMESCALES DURING THE HE-SUBFLASHES BEFORE THE 
BEGINING OF THE EHB. Tduration INDICATES THE LENGTH OF THE 
UNSTABLE STAGES DURING EACH HE-SUBFLASH, AS SHOWN IN FiG. 
[21 Te IS THE E-FOLDING TIME OBTAINED FROM THE NONADIABATIC 
PULSATION COMPUTATIONS. Tgio IS A GLOBAL ESTIMATION OF THE 
CONVECTIVE TURNOVER TIME (SEE TEXT). 



short period p-mode pulsations in LS IV-14°116. In the 
light of our results, LS IV-14°116 would be the first star 
in which pulsations driven by the e-mechanism have been 
detected, thus being the first proof that the e-mechanism 
can indeed drive pulsations in stars. 

2. SCENARIO AND NUMERICAL DETAILS 

To explore the effects of the e-mechanism in a con- 
tracting pre He-core burning subdwarf star, we have com- 
puted the evolution of an initially 1.03 M (n, Z = 0.02 star 
model (see Miller B ertolami et al.ir2008l for details about 
this sequence) that departs from the red giant branch 
just before the He-core flash, and experiences the flash 
as a shallow mixing hot flash er of 0.47378 M(t, (mild He- 
enrichment, Y = 0.5978; see iLanz et al.l [20041 for a def- 
inition). Stellar structure and evolution equations have 
been solved with LPCDDE stellar evolution code as in our 
previous works (Miller Bertolami et al. 2008 and refer- 
ences therein). Although our computations have been 
performed within the hot-flasher scenario, the existence 
of He-shell subflashes in both the hot-flasher and merger 
scenarios makes our approach valid regardless the actual 
origin of LS IV-14°116. 

As seen in Fig. [l] after the primary He-core flash (not 
shown in the figures) the star contracts towards the He- 
core burning phase through a series of He-shell subflashes 
(loops in Fig. [T]), evolving through the region in the 
TeS — g diagram occupied by LS IV-14°116. For this se- 
quence we analyzed the pulsational stability properties 
of more than 1000 stellar structure models, one every 
five stellar evolution timesteps, covering the whole evo- 
lution from the primary He-core flash to the final set- 
tling on the He-core burning phase. For each model we 
restricted ourselves to explore the £ = 1 g-modes in the 
period range from 500 to 6000 s. The stability anal- 
ysis of the stellar models was carried out by means of 
the linear , nonradial, r i onadia batic pulsation code de- 
scribed in ICorsico et al.l ()2006f ) with the inclusion of the 
e-mechanism mode driving opera t ing in He-burning lay- 
ers as described in ICorsico et al.l ()2009() . Two standard 
simplifications have been introduced. We have adopted 
the "frozen-in convection" approximation in the stabil- 
ity analysis and we have assumed that dS/dt = in the 
non-perturbed background model adopted in the compu- 
tations. While the latter is not rigorously true in con- 
tracting stellar models, the timescale of the He-subflashes 
(~ 1000 — 10000 yr) is nonetheless much longer than the 
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Fig. 1. — Evolution of our sequence before settling on the He-core 
burning phase. Thick line indicates pre-EHB evolution while thin 
line corresponds to the evolution after the EHB. Colored points 
indicate the location of models for which unstable modes have been 
found. These unstable modes are caused by the e-mechanism acting 
in the unstable He-burning subflashes that follow the primary He- 
core flash. Color coding indicates the number of unstable modes 
found in eac h model. Boxe s indi cate the locati on of LS IV-14°116 
according to INaslim et al.l l|20l0l ') (cyan) and ijGreen et al.l 120111 ) 
(grey). 

oscilation periods (500 — 6000 s), and this approximation 
is correct. To assess the validity of "frozen-in convection" 
approximation we computed the turn-over timescale of 
convection at the peak of the helium flashes as 

7-gio = / , (1) 

where i?basc and i?top are the radial coordinates of the 
boundaries of the convective zone and wmlt is the veloc- 
ity of the convective elements within the mixing length 
theory of convection. As seen in the fourth column of 
Table [l] the convective turnover times are much longer 
than the unstable periods obtained in this work (see Fig. 
[2). Thus, the assumption of "frozen-in convection" is 
largely justified. 

3. RESULTS 

In Fig. [T]we can see how the star becames recurrently 
a pulsator as it evolves towards the EHB. Our sequence 
of models is very close to the location of LS IV-14°116 
in the T^g — g diagram. Note that no attempt to fit the 
exact location has been done, as we have just used a se- 
quence from our previous work ([Miller Bertolami et al.l 
[2008). It would be easy to attain a better agreement by 
considering a sequenc e with a lower He-abundance o r ear- 
lier He-core flash f see [Miller Bertolami et al.ll2008l ). Fig. 
[21 shows the evolution of the pulsation periods of our He- 
sdB sequence of models. In the upper panel we depict the 
situation found for the complete pre-EHB phase. Note 
that, each time a subflash develops, we obtain unstable 
g-modes, with the periods getting shorter as the star ex- 
periences subsequent, deeper, subflashes. The bottom 
panel illustrates the situation for the first and strongest 
flash of the series and shows the derived e-folding times 
for each mode. Notably, modes are continuously desta- 
bilized during the whole flash, with e-folding times (te) 
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Fig. 2. — Upper Panel: Temporal evolution of the pulsation pe- 
riods of our He-sdB models during the He-subflashes before the 
EHB. Thick points indicate those periods found to be formally un- 
stable in our computations. Lower Panel: The situation for the 
first and most intense subflash. Colored points indicate unstable 
modes with color coding indicating the e-folding time of each un- 
stable mode. 



of the excited modes getting shorter than the duration of 
the pulse (rduration): indicating that these modes should 
attain observable amplitudes (see also Table [1]) . 

While the whole evolution from the primary He-flash 
to the EHB takes about 2 Myr, the unstable stages add 
up to ^ 68000 yr (see Table [Ij . It is worth noting that 
this refers to the formally unstable stage, including some 
stages with somewhat long e-folding times, and that the 
actual unstable phases might be shorter. Then, within 
this picture, we would expect the variable to non variable 
ratio of He-sdB to be lower than 1/30. It is worth noting 
that this is consistent with the fact that only one He- 
sdB star, o ut of 23 studied objects, has been found to be 
a variable (lAhmad & Jcffery 2005; Ahmad ct al. 2004). 
As shown in Fig. [21 most unstable modes exhibit periods 
in the range 600 — 2000 s close to the periods observed in 
LS IV-14°116 (1954-5084 s), and definitely much longer 
than the periods displayed by sdBVr stars in the same 
region of the Tcff — g plane. In particular, it is worth not- 
ing that our models are able to reproduce the dominant 
periodicity observed in LS IV-14°116 (P = 1954 s). 

Our computations show that not all modes with pe- 
riods within the unstable range are excited. In fact, 
only modes with significant amplitudes in the He-burning 
shell are destabilized. This can be seen in Fig. [3] where 
the real part of the eigenfunction 6T/T and the differen- 
tial work {dW/ dr) for two modes with consecutive radial 
orders {k = 26; stable and k = 25; unstable) are shown. 
This behavior, in which the e-mechanism acts as a narrow 
band filter exciting only modes with large amplitudes in 
the narrow He-burning shell region, has been reported in 
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Fig. 3. — Lower Panel: Values of 5T/T for two consecutive 
modes during the first He-subfiash, fc = 25 (unstable) and fc = 26 
(stable). For the unstable mode, the large amplitude of ST/T 
in the region of the He-burning shell can be easily appreciated. 
Upper Panel: Differential work (dW/dr) for the modes fc = 25 and 
fc = 26. Note the strong driving contrilsution of the burning shell 
{dW/dr > 0) in the case of the unstable mode. In the interest of 
comparison we have arbitrarily magnified 5T/T and dW/dr for the 
stable mode (dotted line). Brown curls indicate the location of the 
He-flash driven convective zone. 



previous works (jKawaler et al.11198 6': 'Cors ico et al.l2009l) . 

It has been mentioned (jGreen~e t al. 20111) that the ef- 
fect of the enhancement of He in nonadiabatic compu- 
tations in the context of the K-mechanism in sdB stars 
might be related with the pulsations observed in LS IV- 
14°116. As a consequence, once might wonder to which 
extent our results are related with the /c-mechanism. 
Also, one might wonder to which extent the so called 
"convective blocking" mechanism (Pesncll 1987) at the 
He-shell flash convective zone might help to destabilize 
the modes. In order to check our results, we have car- 
ried out explicit tests by switching off the derivatives 
of ehc in our nonadiabatic computations. As expected, 
when the e-mechanism is turned off, all unstable modes 
shown in Fig. [51 become stable, indicating that the driv- 
ing mechanism behind the unstable modes is indeed the 
e-mechanism alone. 

As mentioned in previous paragraphs, our results sug- 
gest that the periods of unstable modes become shorter 
when the star experiences deeper subflashes, that is, 
when the burning shells move towards the core. One 
might then wonder where should the shell be located in 
order to excite modes with periods embracing the whole 
range of periods observed in LS IV-14°116. In this con- 
nection, we performed a numerical experiment by artifi- 
cially shifting outwards the location of the burning shell 
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Fig. 4. — Effect of shifting the location of the burning shell on the 
periods of unstable modes in the first subflash. Blue points indicate 
modes destabilized when the burning shell is artificially shifted to 
logg ~ —0.37. Green points show the unstable modes in our stan- 
dard sequence (burning shell at logtj ^ —0.21; see Fig. (Sjl. Lines 
indicate the periods detected in LS IV-14°116, with their width be - 
ing propotional to the amplitude of each mode 1 Green et al.|[20TlI) . 

in our nonadiabatic computations. As shown in Fig. 21 
it would be enough to shift the location of the shell to 
logg =~ —0.37 {q :— l — m{r)/Mi,) to excite modes with 
periods covering the complete range of observed peri- 
ods of LS IV-14°116. Interestingly enough, this location 
is completely within the range of locations of the burn- 
in g shells in the merger scenario (see figures 2 and 3 
of iSaio fc Jeffervl l2000). Whether this means that LS 
IV-14°116 is a post-merger object or that He-shell sub- 
flashes develop differently from that current stellar evo- 
lution models predict (as suggested by recent 3D hydro- 
simulations of the event; iMocak et al.l 120091 ) should be 
object of further study. 

4. CONCLUSIONS 

We have shown that the e-mechanism is able to excite 
pulsations in the He-sdB star LS IV-14°116. Although 
our computations have been performed within the hot- 
flasher scenario, the existence of He-shell subflashes in 



both the hot-flasher and merger scenarios makes our 
conclusions valid regardless the actual origin of LS IV- 
14°116. 

Our results strongly suggest that long period g-mode 
pulsations exhibited by LS IV-14°116 are driven by the 
e-mechanism. If so, it is the first confirmation that the 
e-mechanism can drive pulsations in stars. Interestingly, 
if LS IV-14°116 is a star in a fast evolutionary phase 
then iron and nickel might not have had enough time 
to accumulate in the driving region and drive p-mode 
pulsations, thus explaining also the intriguing absence of 
these modes in LS IV-14°116 lightcurve. 

The fact that our model covers only partially the range 
of periods observed in LS IV-14°116 might be related to 
the possibility that the star is a post-merger object, in 
which the locations of the He-flash shells are different 
from those of hot-fiasher objects. In fact, the numeri- 
cal experiment performed in the previous section suggest 
that nonadiabatic computations of the e-mechanism on 
post-merger stellar models might show better agreement 
with the observed periods of LS IV-14°116. Computa- 
tions of the e-mechanism in post-merger He-sdB stellar 
models would be very valuable. 

In short, we have shown that the LS IV-14°116 atmo- 
spheric values (Teff, g and nno), the existence of long pe- 
riod g-moAe pulsations, and possibly the absence of the 
classic K-mechanism p-mode pulsations can be simulta- 
neously explained if LS IV- 14° 116 is as star undergoing 
a He-shell fiash in its way towards the EHB. 
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